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Abstract We present a study of the inclusive photon spectra from 5.9 million J /ip decays collected with the 
KEDR detector at the VEPP-4M e + e~ collider. We measure the branching fraction of radiative decay J/tp^- 
t?c7, Vc width and mass. Our preliminary results are: M(rj c ) = 2979.4±1.5±1.9 MeV/c 2 , T(?7 C ) = 27.8±5.1±3.3 
MeV, B{J/iP^ricl) = (2.34±0.15±0.40)%. 
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1 Introduction 



J/ip — > ?7 C 7 decay is an Ml radiative transition 
with photon energy about 114 MeV, and a relatively 
large branching fraction about 2%. This is a transi- 
tion between IS states of charmonium, which rate can 
be easily calculated in potential models. However, for 
a long time a significant difference between theoreti- 
cal predictions and experimental data was observed. 
Until the last year the PDG value of this branching 
was based on a single measurement [l| of Crystal Ball 
using an inclusive spectrum of photons in the J / ip de- 



cays with the result B{J/i> -> ^7) = (1.27 ±0.36)%, 
approximately twice less than the theoretical predic- 
tions (2.9±O.5)%0, (2.O±O.8)%0, (2.4±O.3)%0, 
(1.6±1.0)%@, (2.1±O.4)%0, (3.05±0.07)%@. For 
comparison, in El transitions of charmonium the cor- 
responding difference is (20-30)% only!?!], although 
the matrix element has more complex structure. 

In 2009 CLEO published^ a new measurement 
of this branching with the value B(J/ip — > ^7) = 
(1.98±0.09±0.30)%, which is closer to theory. CLEO 
analysis used exclusive rj c decays and found that the 
line shape of this transition is asymmetric. A Breit- 
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Wigner function (BW) alone provides a poor fit to 
data. A Breit- Wigner modified by Et (according 
to theoretical predictions for electromagnetic dipole 
transitions) improves the fit near a resonance, but 
gives a diverging tail at higher photon energies. To 
suppress this tail CLEO used an additional exponen- 
tial factor. 



2 Photon spectrum 

The photon spectrum in J/ip — > r] c "f decay is given 
by the formula 



du) 



- Nj,B 



*W 2 >t ffV) e{u')g{uM) m 



where B is the decay branching fraction, g(oj,oj') is a 
calorimeter response function, e(oj) is a photon detec- 
tion efficiency, dT(oj)/doj is a true photon spectrum, 

and r r?c7 = dojdT(oj) / dio is the full decay width. The 

true photon spectrum can be written as Q 



dTjoj) 
du 



^a^-oj 3 \M\ 2 BW(oj). (2) 
3 mi 



Here M —< rj c \j (u>r/2)\J/ip > is the matrix element 
of the transition, ja(x) = sin(x)/x, e c and m c are 
c-quark charge and mass. 

We see that if oj tends to zero, the matrix ele- 
ment tends to unity. CLEO used in their fit \M\ 2 = 
exp(—u> 2 /8/3 2 ) with j3 = 65 MeV, but such a form of 



matrix element squared is valid for harmonic oscilla- 
tor wave functions only. In all other potentials \M\ 2 
dependence will be proportional to oj~ n . 

Therefore we tried to fit the CLEO data using 
another line shape: at photon energy oj near reso- 
nance the decay probability dT/duj is proportional 
to oj 3 , but at higher energies the factor oj 3 is re- 
placed with oj. We found that the function dT/doj ~ 
oj 3 oj 2 /(ojoj + (oj — oj ) 2 ))BW(oj), where oj is the most 
probable transition energy, is also suitable. In Fig. [TJ 
our fit of CLEO data using such a function is shown. 
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Fig. 1. Fit of CLEO data using dT/doj ~ 
oj 3 ojq / {ojoio + (oj — o)o) 2 )BW(oj). 



In Table [T] results of fits with the CLEO function and 
our function are shown. We see that confidence levels 
are good in both cases. 



Table 1. Fit results of CLEO data using various decay probability functions. 



dV/dw 


M(r)c),MeV/c 2 


r( Vc ),MeV 


N EXC 

"is 


X 2 /NDF(C.L.) 


~ id 3 exp(-uj 2 /8l3 2 )BW(ij) 
~ uj 3 uj 2 / (ujujo + (uj — ujo) 2 )BW (uj) 


2982.4±0.7 
2981.8±0.5 


32.5 ±1.8 
33.6±1.9 


6142 ±430 
6494 ±362 


38.0/38 (0.47) 
39.1/39 (0.47) 



3 KEDR data 



The experiment was performed at the KEDR 
detector of the VEPP-4M collider 0. In Fig. [2] 
a layout of the VEPP-4M collider is shown. It op- 
erates at a peak luminosity about 1.5- 10 30 cm~ 2 s _1 
near the J/ip peak energy. Two methods of beam en- 
ergy determination are used: resonant depolarization 
with accuracy of 8 ± 30 keV and IR-light Compton 
backscattering with accuracy ~ 100 keV 

The view of the KEDR detector is shown in Fig. [3] 
Subsystems are listed in the figure. In this analy- 
sis a barrel liquid krypton calorimeter, drift chamber 
(DC) and time-of-flight (ToF) scintillator counters 



were of decisive importance. Luminosity was mea- 
sured by single Bremsstrahlung online and by small 
angle Bhabha scattering offline. 

This analysis is based on a data sample of 1.52 ± 
0.08 pb" 1 collected at the J /ip peak. Three J/ip scans 
were performed; a measured beam energy spread was 
used for the calculation of the number of J /ip pro- 
duced: ^ = (5.9±0.3)-10 6 . 

Event selection was done in two steps. At the first 
step multihadron decays of J/ijj were selected. To do 
this the following cuts were applied: total energy in 
clusters in calorimeters is greater than 0.8 GeV; at 
least 4 clusters with energy greater than 30 MeV in 
calorimeters are detected; at least one central track 
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in the drift chamber is reconstructed; there are no 
muon tubes activated in the third layer. These cuts 
effectively suppress background from the cosmic rays, 
beam-gas interactions and Bhabha events. 




Fig. 2. VEPP-4M collider. 




Fig. 3. KEDR detector. 1-vacuum cham- 
ber, 2-vertex detector, 3-drift chamber, 4- 
threshold aerogel counters, 5-ToF-counters, 6- 
liquid krypton calorimeter, 7-superconducting 
coil (0.65 T), 8-magnet yoke, 9-muon tubes, 
10-CsI-calorimeter, 11-compensation solenoid, 
12-VEPP-4M quadrupole. 

At the second step photons in these events were 
identified. A photon is a cluster in the liquid krypton 
calorimeter without reconstructed tracks in the drift 
chamber attached to it and without ToF scintillator 
counters activated before the cluster. According to 
the simulation, the photon detection efficiency with 
the above mentioned cuts is about 34%. 



4 Data analysis 

The inclusive spectra of clusters are shown in 
Fig. Large peaks in the spectra correspond to 
minimal ionizing particles crossing the calorimeter. 
The spectrum of charged particles was used to re- 
ject charged particles which were detected as neutrals 
from the photon spectrum during fitting procedure. 



xlO 3 




co, MeV 

Fig. 4. Inclusive spectra in the liquid kryp- 
ton calorimeter: open squares - all clusters; 
shaded circles - charged particles, when the 
ToF scintillator counter is activated before a 
cluster and a DC track is attached to cluster; 
open circles - photons. 

In Fig. [5] a fit to our data is shown. 
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Fig. 5. Fit of the inclusive photon spectrum 
with f(co) = loq/(ujloq + (uj — cjq) 2 )- 



The inclusive photon spectrum was fit by the sum of 
the signal having a shape dT/dui ~ u 3 f(uj)BW(uj), 
convolved with the calorimeter response function 
(Novosibirsk PDF - asymmetric Gaussian with <r E = 
8.8 MeV at 110 MeV and a=-0.3), and background. 
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Fig. 6. Resulting signal line shapes for various 
dT/du functions. 

Background was taken in the form of exponent of a 
second-order polynomial plus a scaled histogram of 



charged particles. A relativistic Breit-Wigner func- 
tion was used in the fit. 

We also tried to fit our data with various line 
shapes. We used BW alone, BW-w 3 : CLEO function, 
and our function. The resulting signal line shapes for 
these functions are shown in Fig. [6] The BW func- 
tion alone gives a shifted value of r) c mass compared 
to other functions: M(r? c ) = 2975.9 ±1.3 MeV/c 2 . 
BW-uj 3 gives a better form near the resonance, but 
leads to the unphysical diverging tail at higher pho- 
ton energies. The last two functions give close results 
of fits, and their difference can be used for estimation 
of a systematic error appearing due to the unknown 
line shape. In Table [5] results of fits of KEDR data 
with the CLEO function and our function are listed. 
We see that confidence levels are good in both cases. 



Table 2. Fit results of KEDR data using various decay probability functions. 

dT/duj M(r) c ),MeV/c 2 r( Vc ),MeV B(J/ip) -+7t?c), % X 2 /NDF{C.L.) 

~uj 3 exp(-uj 2 /8l3 2 )BW(uj), p = 65MeV 2979. 7±1. 6 26.9±4.8 2.19±0.28 84.1/78 (0.30) 

~aj 3 wJ/(wuo + (w-i*) 2 )BfF(u) 2979. 4±1. 5 27.8±5.1 2.36±0.34 83.9/78 (0.30) 



5 Systematic errors 

In Table [3] systematic errors of our measurements 
are listed. Systematic error appearing due to the un- 
known line shape was estimated using various line 
shapes and fixing r) c width to the current PDG val- 
ues (27.4±2.9 MeV). Systematic error related to the 
background subtraction was estimated by the taking 



a polynomial of the third order instead of the second 
in the fit, and varying ranges of the fit. Luminosity 
measurement error was evaluated from the difference 
between the measured and calculated number of the 
multihadron J/tp decays. Photon detection efficiency 
is conservatively estimated as 10%. Photon energy 
scale calibration was made comparing a 7r° peak po- 
sition in simulation and data. 



Table 3. Systematic errors. 



Systematic error 



M(r! c ),MeV/c 2 



r( Vc ),MeV 



e(J/^)^ 7 r;c),% 



Line shape 

Background subtraction 
Luminosity measurement 
Photon efficiency 
Photon energy scale 



1.5 
0.5 



1.1 



2.5 
2.2 



0.25 
0.18 
0.11 
0.23 



Total 



1.9 



3.3 



0.40 



6 Results and conclusions 

A new direct measurement of J/ip — > rj c ^ decay 
was performed. We measured rj c mass, width, and 
branching fraction of J/ip — > r/ c 7 decay. Our prelimi- 
nary results are: 

M(j7c) = 2979.4±1.5±1.9 MeV/c 2 , 

r(?7 c ) = 27.8±5.1±3.3 MeV, 

B(J/V> -> f?c7) = (2.34 ± 0.15 ± 0.40)%. 



In Fig. [7] results of our measurements, compared 
with the other experiments and the current PDG val- 
ues, are shown. The branching fraction of this decay 
is sensitive to the line shape, which is asymmetric, 
and it should be taken into account during a mea- 
surement. 

We see that our results on rj c mass and width are 
consistent with the current PDG values. Our result 
for B{J/ip — > 77 c 7) is consistent with that of CLEO, 
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Fig. 7. 77c mass, width and B(J/i/j —> ^7) measurements. The theoretical predictions for the branching 
fraction are also shown (squares). A small error of Eichten 08 prediction is a consequence of the fact that the 
authors take into account the experimental uncertainty of the J/ip total width only and ignore all theoretical 
uncertainties. 
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